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This  paper  deals  with  the  development  and  analysis  of  well-posed  models  and  compu¬ 
tational  algorithms  for  control  of  a  class  of  partial  differential  equations  that  describe  the 
motions  of  thermo- viscoelastic  structures.  We  first  present  an  abstract  “state  space”  frame¬ 
work  and  a  general  well-posedness  result  that  can  be  applied  to  a  large  class  of  thermo-elastic 
and  thermo-viscoelastic  models.  This  state  space  framework  is  used  in  the  development  of 
a  computational  scheme  to  be  used  in  the  solution  of  an  LQR  control  problem.  A  detailed 
convergence  proof  is  provided  for  the  viscoelastic  model  and  several  numerical  results  are 
presented  to  illustrate  the  theory  and  to  analyze  problems  for  which  the  theory  is  incomplete. 
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1.  Introduction.  During  the  past  few  years  considerable  attention  has  been  de¬ 
voted  to  the  development  of  smart  materials  and  structures  (see  [B] ) .  One  approach 
to  this  class  of  problems  is  to  use  shape  memory  alloys  as  actuators  in  active  control 
designs.  These  alloys  are  best  described  by  thermo-mechanical  models  consisting  of 
coupled  (and  nonlinear)  hyperbolic  and  parabolic  partial  differential  equations.  The 
development  of  computational  algorithms  for  designing  controllers  for  such  systems  is 
an  immensely  complex  problem  and  the  subject  of  several  ongoing  research  projects. 
In  addition  to  the  obvious  difficulties  related  to  the  nonlinearities,  the  basic  thermo¬ 
elastic  coupling  often  gives  rise  to  nonstandard  mathematical  models  and  leads  to 
several  problems  in  developing  computational  algorithms  for  control.  Therefore,  the 
computational  methods  for  controlling  a  linear  thermo-elastic  system  may  be  viewed 
as  a  first  step  toward  the  ultimate  nonlinear  problem.  With  this  motivation  in  mind, 
we  consider  the  problem  of  controlling  a  class  of  coupled  partial  differential  equations 
that  describe  the  linearized  motions  of  a  thermo-mechanical  structure.  The  basic 
approach  is  to  combine  approximation  theory  with  state  space  modelling  to  develop 
convergent  computational  algorithms  for  LQR  control  designs. 

In  this  paper  we  consider  the  questions  of  well- posed  ness  and  convergence  of  ap¬ 
proximation  schemes  for  a  class  of  abstract  linear  systems  of  the  form 

z(t)  =  Az(t)  +  Bu(i).  z(  0)  =  r„  (1-1) 

on  a  Hilbert  space  Z.  The  main  concern  of  this  paper  is  with  a  general  class  of  partial 
functional  differential  equations  (PFDEs)  arising  in  the  modeling  of  viscoelastic  and 


thermo-viscoelastic  systems,  for  example  (see  [BCLM],  [BMC],  [MH]),  coupled 
equations  of  the  form 

=  I  {T^'-r)  +  Sj(s)h(t  +  s-z)i\ 

Q 

~  1~dx°^i'X)  +  (1-2) 

±0it.r}  =  ,±,0il.z)--,0o^gly{,.T)  (1.3) 

where  y  represents  displacement  and  0  is  the  deviation  from  the  reference  temperature 
Oq.  Equations  of  viscoelasticity  (e.g..  (1.2)  with  7  =  0)  have  a  special  structure  which 
has  been  used  by  Fabiano  and  Ito  [FI]  to  formulate  a  general  well-posedness  theorem 
and  convergence  results.  Observe  that  (1.2)  with  7  =  0  can  be  written  as 

d' 2  02  f°  1 

Cr^y^"r)  =  d72  T^  +  J  My^  +  ^x'fds  +b{x)u(t ),  (1.4) 

or.  in  abstract  form,  as 

y  +  A  ry  +  J  g(s)y{t  +  s)ds  =/(/),  (1.5) 

where  A  is  a  positive  definite,  self-adjoint,  closed  linear  operator  on  a  Hilbert  space 
V  .  In  [FI]  Fabiano  and  Ito  consider  equations  of  this  form  with  singular  kernels 
(i.e.,  g  €  Ll(  —  r.O))  and  establish  well-posedness  when  the  state  space  is  taken  to  be 
P(.4I/2)  x  Y  x  L2g{  —  r,  0;  V{AX^2)).  In  this  paper,  we  also  consider  equations  of  the 
form  (1.5).  but  the  approach  we  take  also  applies  to  the  thermo- viscoelastic  equations 
(1.2)  -  (1.3)  which  cannot  be  written  in  the  form  of  (1.5). 

In  Section  2  we  develop  an  abstract  framework  and  a  generalized  well-posedness 
theorem  which  we  apply  to  the  thermo- viscoelastic  system  (1.2)  -  (1.3)  with  zero 
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boundary  conditions  and  to  the  general  viscoelastic  system  given  by  equation  (1.5). 
Our  approach  allows  a  singular  kernel,  and  it  also  has  the  advantage  that  it  does  not 
require  explicit  knowledge  of  the  domain  of  .41//2  in  order  to  write  down  the  state 
space.  This  property  can  be  useful  in  applications  where  A1^2  is  not  a  differential 
operator.  We  also  remark  that  our  general  framework  can  be  applied  to  certain  finite 
delay  systems  similar  in  form  to  the  infinite  delay  systems  considered  by  Miller  and 
Desch  in  [MD],  Miller  and  Desch  prove  well-posedness  for  a  class  of  equations  in 
which  the  kernel  is  completely  monotonic. 

Approximation  of  such  systems  generally  consists  of  two  steps:  first  approximate 
the  spatial  variable  (e.g..  by  means  of  splines)  to  reduce  the  system  to  a  hereditary 
differential  system  on  Rn,  then  approximate  the  "history"  or  "memory"  term  (i.e.. 
the  integral  term  in  (1.5)).  In  this  paper  we  will  use  a  variation  introduced  by  Fabiano 
and  Ito  ([FI])  of  the  averaging  scheme  considered  by  Banks  and  Burns  ([BB])  for 
the  second  stage.  The  idea  of  the  "AVE"  scheme  is  essentially  to  approximate  the 
kernel  g[s)  bv  a  step  function:  partition  [ — r.  0]  into  M  subintervals  and  take  the 
integral  average  in  each  subinterval.  Fabiario  and  Ito  show  that  the  approximation 
scheme  converges  for  an  V  kernel  using  a  uniform  partition  of  [  —  r.  0],  but  they  give 
numerical  results  which  indicate  that  a  different  partition  using  a  finer  mesh  near  the 
singularity  at  zero  yields  much  faster  convergence.  In  Section  3  we  modify  the  proof 
given  by  Fabiano  and  Ito  for  singular  kernels  and  a  uniform  mesh  to  include  singular 
kernels  and  the  non-uniform  mesh. 
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Although  we  prove  convergence  in  this  paper  only  for  the  abstract  viscoelastic 
system,  the  proof  we  give  can  be  modified  to  include  the  thermo- viscoelastic  model. 
In  Section  4  we  give  some  numerical  results  comparing  the  viscoelastic  and  thermo¬ 
viscoelastic  models. 

We  will  use  the  following  notation.  For  a  function  g  £  Ll(a,b\Z ),  we  denote 
by  L2g(a,b:Z)  the  set  j/  £  L2[a.b\Z)  g{s)  ||/(^)]||  ds  <  ocj.  We  denote  by  the 
symbol  H[(a.  6;  Z )  the  set.  of  all  //'  functions  which  vanish  at  the  left  end-point  of 
the  interval;  i.e.,  lPL{a.b :  Z)  =  {/  £  //’(«,  b:  Z)  |  /(«)  =  0}.  Sin  ilarly,  6;  Z)  = 

{/  £  Hl{a.b:Z)  \  f{b)  -  0}.  For  a  function  x  :  [-r.a)  — +  A',  r.  a  >  0.  the  symbol  xt 
for  t  £  [O.a)  represents  the  function  xt  :  [  — r,0]  — >  A’  defined  by  xt(s)  =  x(t  +  s).  If  .4 
is  the  infinitesimal  generator  of  a  C0  semigroup  l'(-)  on  a  Hilbert  space  Z  satisfying 
\\T(t) \\z  <  M e3t .  then  we  write  .4  £  G{M,3).  Finally.  zn-^>z  means  that  zn 
converges  strongly  to  r. 

2.  Well-Posedness.  A  standard  technique  for  establishing  well-posedness  of  a 
system  governed  by  a  PFDE  is  to  cast  the  problem  in  the  form  of  (1.1)  and  show 
that  A  generates  a  Co  semigroup  on  2,  for  example,  by  means  of  the  Lumer-Phillips 
theorem  (see  [P] ) .  We  will  use  the  following  version  of  this  theorem: 

T II EO  Ft  EM  2.1.  Let  .4  be  a  closed  densely  defined  linear  operator  on  a  Hilbert  space 
H.  If  there  exists  3  £  R  such  that  (Ax,x)  <  3(x,x)  for  all  x  £  'D(A).  and  7v(A 07  —  .4) 
is  dense  in  H  for  some  A0  >  3.  then  .4  is  the  infinitesimal  generator  of  a  Co  semigroup 
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T(t)  on  H  satisfying  ||7'(<)||  <  edt. 


2-1.  A  General  Theorem  on  Well-Posedness.  Suppose  that  A\  Y,  0  and  IV 
are  Hilbert  spaces,  and  set  Z@  =  X  x  Y  x  0  x  IV.  Let  5  be  a  subspace  of  Y.  and 
suppose  we  have  the  following  linear  operators: 


Ao 

■  D(.40)  C  Y  -  V, 

.4,  :  I >04,)  C  A'  -  V. 

G\ 

:  D(G]  )C0-»  Y. 

Gi  :  V(G2)  C  Y  —  0. 

63 

■  X>(6'3)  C0-0. 

C\  ■  V(C'i)  C  IV  -  Y. 

D  : 

V{D)  C  IV  -»  IV. 

i  ■  -V  ->  U  '.  j  :  5  - 

.V. 


Define  .4,  C  and  G  by  A  —  .4q.4j.  C  —  .4qCi  and  G  =  A^G^.  and  define  Fq  by 


T>(F@)  = 


Fq  9 


0  1  €  .V  x  0  x  IV 
vc 


X  e  p( .40,0  6  V(G,).w  e  V(C\ 
AjX  +  G\9  +  Ci w  G  V(A0) 


=  .40(.4i.r  +  GiO  +  Cite). 


Define  Aq  bv 


f  x\ 


V(Ae)  = 


y 

e 


€  Zq 


Aq 


(  x  \ 

y 

o 

\wj 


\WJ 
( 


y  G  5  n  V(G2).9  g  V(G3). 

weV(D)Je)  ex>(F©) 


jy  \ 

X 

Fq  [  9 

w 

02y  +  G30 
\  ijy  +  Dw  ) 


Finally,  suppose  that  j  is  injective  and,  for  A  G  p{D)  D  p(G3).  define  L\  :  V(LX)  C 
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A'  -  V  by 


V{Lx)  =  {x  €  K(j)  I  {T,\j~lx,{\ I  -  G3)-lG2\]-'x,(\I  -  D)~'i\x)T  G  V(Ae)}  . 

L\x  =  A 2j~xx  —  F&  (  (A/  -  G3)_1  G2Aj  'x  |  . 

V  (A/  -  D)~li\x  ) 

We  are  now  ready  to  state  the  main  result  of  this  chapter. 


Theorem  2.2.  Suppose 

(1)  ?  and  j~]  are  continuous. 

(2)  T>( .4 )  is  dense  in  A.  5  C  U[G?)  and  S  is  dense  in  Y.  T>[G)C\  V{G3)  is  dense  in 
0.  and  D,C)  1-1  F(D)  is  dense  in  IT. 

(3)  j(S)  is  closed  in  A. 

(4)  Fq.  6' 3  and  D  are  closed. 

(5)  for  y  G  S.  \\G2y\\e  <  k  ||ji/||A.  for  some  k  >  0. 

(6)  there  exists  3  G  R  such  that  (A@z.z)z&  <  3(z.z)z&  for  all  z  G  P(^40). 

(7)  there  exists  A0  >  d.  A0  G  p{D)  fl  p{G3),  such  that  TZ{L\0)  is  dense  in  Y.  and 

(8)  (A 0I  -  D)[D{C)nV{D)}  is  dense  in  IV’,  and  {\0I  ~  G3)[V{G)  DV(G3)}  is  dense 
in  0. 

Then  Aq  is  the  infinitesimal  generator  of  a  C0  semigroup  So(t)  on  Z©  satisfying 
\\Se(t)\\<e^. 


PROOF:  Set  V  =  r>(v4)x5x(P(G)nr>(G3))x(D(C)nD(£>)).  Then  V  C  P^ejand 

(  Xn  \ 

V  is  dense  in  Z©.  so  T)(Aq  )  is  dense  in  Z®.  For  n  =  1,2 . let 


V  n 

On 

\  Wn  J 


G  V(A&).  and 
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/  Xn  \ 

^  I  \ 

fx") 

(A 

suppose 

Vn  1 
0n  , 

- 

y 

e 

and  Aq 

!  Vn 

1  0"  , 

= 

Vn 

In 

V! 

1  T 

\w*) 

\wj 

l  U’n  J 

WJ 

w 

as  n  — *  oc.  Then 

yn  £  S  and  jyn  =  pn  —>  p.  Since  j(S)  is  closed,  there  exists  y  £  S  such  that  jy  —  p. 
But  j~l  is  bounded,  so  ||y- j/||  <  ||j_1j|  •  ||f  ~  Fnll  +  l|i/«  -  //II  — ♦  0  as  n  — >  oc. 
Therefore,  y  -  y:  i.e.,  y  £  S  and  jy  —  p.  Now,  f  0n  j  £  P(F©).  (  0n  J  — ►  (  0 

f  In  \  .  ,  .  f  X\ 

and  F©  I  0n  \  =  vn  — *  V  as  n  — »  oc.  Since  F©  is  closed.  0  j  £  V(Fq)  and 

f  X\ 

Fq  I  0  =  v.  Since  jyn  — »  jy  and  i  is  continuous,  we  have  ijyn  — ►  ijy.  Wo  also 

have  ijyn  +  Dup  — >  h.  Thus.  Dwn  —>  h  —  ijy.  But  D  is  closed,  so  ic  £  D{D)  and 
Div  =  h  —  ijy  which  implies  that  ijy  -f  Dw  —  h.  Next.  F,  — »  0.  0n  £  T>(G3)  and 


G2yn  +  G30n  =  In  T  By  (5).  ||G2(S/n  ~>j)  II©  <  *  Wji'Jn  ~  2/ )  I !  A'  =  k\\jVn  ~  3V  II.V 
0.  so  G2yn  —>  G2 y  which  implies  that  G30n  — ►  o  —  G2y ■  Since  G3  is  closed.  0  £  V(G3) 
and  G36  =  ~)  —  G2y.  or  7  =  G2y  +  G30.  Therefore.  -4©  is  closed.  Finally,  let 


J_  7v(A0/  —  Ao): 


i.e.. 


w 


(vn  A0x  -  jy)x  +  /  riK  X0y  -  Fq  (  0  j  \ 

\  V  W  J  /  Y 

+  {l-  (A0/  -  G3)0  -  G2y)Q  +  (/;,  (A 0/  —  D)w  —  ijy)w  =  0 


/  x  \ 

( 

X 

\ 

for 

all 

1 

y 

0 

£  'D(Aq).  Let 

x  € 

Wo)- 

Then 

(A  0/- 

A0i  ~ 1  x 

G's  )  _  1  G'2  AqJ 

1 

6 

\wj 

V  (Ao/ 

1 

1 

sL 

0 

0 

) 

We), 

so  ( 

b\  AqJ  ~ 1  x  -  F©  f 

( A0/ 

J 

—  0  3  )  ^  ^7  0  Art  /  ^  X 

))  - 

(V.LXox)y 

=  0 

for 

1 

\  V 

(A 

0/  -/})-’ 

?A0x 

all 

x  € 

V{L 

a0),  which  implies 

0  by  (7) 

.  If  x  = 

:  0.  y  = 

0  and  0  = 

0.  then 

/ 


(MV  -D)w)w 

(  x\ 


X  6  V(A), 


0 

0 


=  0  for  ail  w  £  P(C')  D  'D(D),  and  hence  h  —  0  by  (8).  Now  for 
£  V(Aq).  so  (i,p,  A0x)v  =  0  for  all  .r  6  D(/l).  By  (2)  ibis  implies 


w 

that  y?  =  0.  Finally,  for  0  £  P(G’)  fll^Gd). 


/o\ 


0 

0 


£  V(A&),  so  (7,(Aq/  —  O3)0)e  —  0 


Vo/ 

which  implies  that  7  =  0  by  (8).  Therefore.  7v(A 0I  —  Aq)  is  dense  in  Zq.  and  this 
completes  the  proof.  | 


We  wish  to  apply  this  theorem  to  thermo-viscoelastic  systems  (equations  (1.2)  - 
(1.3))  and  to  the  abstract  viscoelastic  system  (equation  (1.5)).  Before  we  proceed 
with  any  examples,  however,  we  will  make  some  general  comments  about  the  kernel 
function  g  and  the  space  IF. 


HYPOTHESIS  2.3.  The  function  g  satisfies  the  following  conditions. 

(1)  g  £  /. 1  ( —  r,  0 ) . 

(2)  g  <  0  and  g1  <  0  on  [  —  r.  0).  and 

(3)  a  =  r  +  J°rg{s)ds  >  0. 


Set  <7o(.s)  =  —-hg(s).  and  suppose  the  space  X  is  given, 
space  Lg(  —  r,0:X)  with  inner  product  given  by 


(w 


ga(s)  (te1(s),ie2(s)) 


We  will  take  W  to  be  the 


A'  ds- 


Define  the  operator  D  by  T>(D)  =  h/^(—  r.  0:  A  ),  D  =  The  following  lemma  is 
proved  in  [FI].  Since  it  is  the  crucial  step  which  allows  a  singular  kernel,  we  reproduce 
its  proof  here. 
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LEMMA  2.4.  The  operator  D  is  dissipative  in  W. 


PROOF:  For  w  e  T>{D), 


( Dw.w)w  =  J  g*{s)  (J^u'{s),w{s))  ds  =  l 


X 


IIU'(-?)IIa'  df- 


Let  e  >  0  and  consider 


h 


u-(s)|| 


A' 


ds 


=  |k’(-£)lli  “  lk’(-r)llA'  j  9o(S 

<  £&»(-«)  IM-<0!!2y  • 


|ir(s)||2v  ds 


Since  iv(—  e)  =  u'(0)  —  f°(  Dw(s)ds  =  —  J^Dw(s)ds,  by  the  Cauchy-Schwarz  in¬ 
equality 

IM-c)||2v  <  f  ~T~\  f  do{s)  ||Z)?c(s)||2v  ds. 

J -c  9q\s)  J -( 

/0  rO  / \ 

-  =  /  — — — ds  <  e.  Thus,  we  obtain 

■«  9o(s)  J-c  ga[s) 

h  <  -}  j  go{s)  ||/?ti;(s)||2v  ds 

for  all  e  >  0.  Therefore  (T >w,w)w  =  liirpjo  h  <  0.  I 
For  h  €  IF,  if  u>(s)  =  es  e~° h(a)do ,  then  €  T>(D)  and  (/  —  T>)u'  =  ft.  Since  D 
is  densely  defined,  D  generates  a  Co  semigroup  of  contractions  on  IF  by  the  Lumer- 
Phillips  Theorem.  In  particular,  D  is  closed,  and  A  £  p{D)  for  all  A  >  0. 


2.2.  A  Thermo- Viscoelastic  System.  We  consider  the  system  governed  by  equa¬ 
tions  (1.2)  -  (1.3)  with  boundary  conditions  given  by  y(t,  0)  =  y{tA)  =  O(t.Q)  = 
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0(1,  1 )  =  0.  If  we  set  w(t,  s,  x)  =  y(t,  x)  —  y(t  +  s,  x)  and  use  ga  as  defined  above,  then 


we  can  rewrite  equation  (1.2)  as 

~y(t,x)  =  y(t,x)+f  ga(s)-^-u'(t,s,x)ds  -  -0(1,x)  +  -b(x)u(t). 

Ol 1  a  ox  [ox  J_T  ox  a  J  o 

Define  the  spaces  X,  Y  and  0  as  follows: 


A  =  //(](0.  1 )  wit  h 

H 

to 

II 

O 

/  ✓,* 

Jo 

Y  =  L2(0. 1 )  with 

(j/i  •  3/2))-  =  <7 

f  7/17/2- 

J 

0 

0  =  L2(0.  1 )  with 

(<v<y9  =  J- 

/  0i  02 

Jo 

Set  Zq  =  A  x  Y  x  0  x  IF  where  IF  =  L2g(—  r.  0:  A  )  is  as  above,  and  let  5  C  V’  be 
given  by  5  =  H^( 0,  1).  Define  the  following  operators: 

■Z>(/to)  =  //,(0.1).  A0y  =  -y'  6  Y. 

a 

P(Ai)  =  .A,  Alx  =  x'eY, 

T>(Gi)  =  0.  C,tf  =  --0eV, 

Q 

V(G2)  =  H'(OA).  ^  =  -iVe0, 

V(G3)  =  Hl( o,  i )  n  H2( 0,1),  G3e  =  k0"  e  0, 

P(C])  =  IF,  C, a-  =  J  ga(s)w'(s)ds  <E  K, 

[zx)(s)  =  x  6  IF,  y:  7/q  =  5  — >  A'  =  77q  is  the  identity  operator. 

Wit h  the  above  definitions,  we  have 

V(A)  =  III  n  H\  T>(C)  =  L2g(-r,0;V(A)),  V(G)  =  H\ 0.1). 
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and  the  operator  Aq  is  given  by 


T>(Ae) 

(* 


\ 

(x\ 

y 

0 

6  Zq 

1 

\w) 

y 

o 

\w 


y  6  Hl,0  6  HI  OH2. 

*'  -  l0  +  S-r  9a{s)iv'{s)(ls  £  77  1 


/  V  \ 

f  £  (X'  “  l0  +  f-rSMw'(s)<ts) 

-^Ooy’  +  KO" 


y  +  T-s 


) 


We  now  verify  the  conditions  of  Theorem  2.2. 


(1)  Clearly  ?  and  j_1  are  continuous. 

(2)  Clearly  V(A)  is  dense  in  ,Y,  5  C  PfGj)  and  S  is  dense  in  V’.  T>(G)  PI  22(G3)  = 
F>(G3)  is  dense  in  0.  and  V{C)  fl  22(72)  =  //^(-r,0:'D(.4))  is  dense  in  IV. 

(3)  Since  j(S)  =  X .  j{S)  is  closed. 

(4)  We  already  know  that  D  is  closed.  It  is  easy  to  see  that  G3  is  densely  de¬ 
fined  and  dissipative,  and  for  any  G  0  there  exists  0  6  Z>(G3)  such  that 
(7  —  G3)0  —  v?  (see  [K.  p.  147]).  Thus.  G3  is  closed.  Let  f  0n  \  £  V(Fq). 


w 
a  d 
(7  dx 


w 


Xn 

On 


0n  |  — ^ ■  (  0  ,  and  Fq  0n  I  =  b’n  — ♦  A-  Observe  that  F@  6n  = 


IL'n 


^x'n  -  ?-6n  -f  ga(s)w'n(s)dsSJ .  Now,  xn  — *  x  in  A'  implies  that  x'n  — ♦  x' 
in  Y .  #n  — »  0  in  0  implies  that  0n  —>  0  in  Y,  and  wn  — »  u;  in  IV  implies  that 
f°r9o(s)K(s)ds  -»  J°rg0{s)w'(s)ds  in  V.  Thus,  x'n  -  +  f°T  g0{s)u>'n(s)ds 

^x'-ie  +  fl9  0{s)w'(s)ds  in  Y .  Since  ,4o  is  closed.  F&  is  closed. 

(5)  It  is  easy  to  check  that  ||G’2J/j|©  =  — —  ||ji/||^.  Set  k  =  >  0. 
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(6)  For  r  =  (x,y,d,w)T  G  P{A&), 

(A&z,z)=q  f  y'x'  +  a  f  (x'  -  -0  +  [  gQ(s)w'[s)ds\  y 

Jo  Jo  crdx  \  o  J_T  J 

+  Yq  J  (-7^o y'  +  n0")0  +  J  ga(s)  y'w'  +  (Du\  w)w 

=  K  { G30,0)q  t  {Pw,  w)w  <  0 
since  6’3  and  D  are  dissipative. 


(7)  We  will  take  A0  =  1.  Let  x  G  V(L\).  If  we  set  w(s)  =  (1  —  es)ix.  then 
(/  —  D)w  —  i. r.  so  (  /  —  D)~xix  =  (1  —  ts)ix  =  (1  —  e5)x.  Note  that 

L\X  =  x  —  —  —  (x1  H - —  (/  —  G'i)  +  J  <7er(-s)(l  —  ts)xds 

Q  d  (  l2°o,r  n  , 

=  ,r - —  QiJ  +  - — -(/  -  G3  x 

a  ax  \  a 


where  a  \  =  — 
a 


l  r 


Q  -  /°r  $(*)(!  -  £Xs)ds 


a  L 


+  f°rg{s)eXsds 


>  0  for  A  > 


0.  Since  .Y  C  Y ,  we  can  think  of  L\  as  being  defined  on  a  subspace  of  Y. 
Observe  that  V(L\ )  =  j y  G  //q(0,  1 )  |  a}y'  +  -^{I  -  G'3)~V  G  H1 1  is  dense 
in  }  .  Define  the  operators  p  and  T2  as  follows: 


V(P)  =  H\  0,1),  Py  =  y', 
V{T2)  =  //q(0,1),  T2y  =  y'. 


Note  that  T\  and  T2  are  adjoint  to  each  other.  With  this  notation  we  can  write 
Li  as  follows: 


P  =  I -~P 

a 


C1  +  -U-G,)- 1 

Q 


TV 
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Since  Gj  =  G 3,  we  have  [(/  —  G 3)  !]*  =  (/  —  G 3)  T  Thus, 


Q 


l;  =  ;  -  -Tj 


I--T, 

a 


a,l  +  —V  -G3y' 


t r 


Ql/  +  2^2(/_G3)-i|  72  =  /„; 

Q 


that  is,  Lj  is  self-adjoint.  Now,  for  y  6  P(I]), 

(L^y.y)y  =  (y.  y)v-  -  o  J  ~  ^a,y#  +  -^(/  -  GsTV^  y 

=  llyliv-  +  Q  J  y'  +  -  G3)- 2/ 

=  1 1  y  1 1  v'  +  or,  \\jy\\2x  +  -*<(/-  ^rly'<y')Y  >  Ml 


since  I  —  G 3  >  0  implies  that  (/  —  G 3)  1  >  0.  Thus.  Tj  is  one-to-one.  Hence 
by  Theorem  13.11  in  [Ft] ,  TZ(Li)  is  dense  in  Y. 

(S)  Since  V(G3)  C  D(G).  we  have  (/  -  G3)  \V(G)  n  V{G3))  =  (I  -  G3)V{G3)  = 
0  from  (4).  Next,  for  h  €  P(C),  if  we  set  w{s)  =  es  J°  c~° h(cr)da,  then 
w  e  HlR{-r,0;V(A))  =  V{C)  n  V{D)  and  (/  -  D)w  =  h.  Thus,  V(C)  C 
(/  —  D)  [D(C)  fl  V(D)],  and  V(C)  is  dense  in  W . 

Since  (1)  -  (S)  hold.  A@  generates  a  Co  semigroup  on  Z©. 

2.3.  The  Abstract  Viscoelastic  System.  Next  we  wish  to  apply  Theorem  2.2 
to  equation  (1.5).  Since  the  space  0  corresponds  to  the  heat  equation  (1.3),  we 
first  formulate  a  version  of  the  theorem  using  only  the  spaces  .V,  V  and  IT.  Set 
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Z  =  .V  x  }'  x  li’,  and  define  the  operator  F  by 


F(F)  —  j  (  ^  )  €  A  x  If  i  G  T-\ A \ ),  w  G  T)(C  i ),  A\x  4-  C\w  G  ’P(-do)  |  . 
F(w)  =  -do(.4ix  +  Citi’), 


define  A  by 


V(A) 


y  G  S,w  €  F{D),  | 

C)e»(n  /• 


and  define  L\  by 


V(LX)  =  {x  G  1l(j)  |  (x.\j~'x,(\I  -  DyUXxf  G  V(A)}  . 
Lxx  =  -  F  (^(A/  _  a*)  ■ 


We  have  the  following  special  case  of  Theorem  2.2. 

THEOREM  2.5.  Assume  that  all  the  hypotheses  of  Theorem  2.2  except  those  involving 
the  space  6  jiold.  In  addition,  assume  that  F  is  closed,  that  A  satisfies  the  inequality 
in  (6),  and  that  there  exists  A0  >  T,  A0  G  p{D),  such  that  7 Z(L\0)  is  dense  in  Y .  Then 
A  generates  a  Co  semigroup  S{t)  on  Z  satisfying  [|S'(<)[|  <  e'3t . 


PROOF:  Let  0  =  {0},  and  define  G \  =  0,  G2  —  0  (with  ^>(^2)  =  Y)  and  G 3  = 


0.  Define  V(Fq)  = 


0  G.YxQx  W 

w  ) 


and  define  Aq  by  V(Aq)  = 


G  Z, 


0 


(J)€D(F)J.Fe(o)  =  f(£). 


y  )eV{A)}  .  A o,~6  =  TAT~l:e 

w 


Z  whose  infinitesimal  generator  is  A.  | 

We  wish  to  apply  Theorem  2.5  to  the  equation  (1.5).  If  we  set  w{$)  =  y  -  ys  and 
substitute  into  (1.5)  we  obtain  an  equation  of  the  form 

y  +  aA  y  +  J  g0[s)u:(s)di.  =  /(/).  (2.1) 

We  assume  that  A  is  a  closed,  densely  defined,  positive,  self-adjoint,  injective  linear 
operator.  As  we  shall  see  in  the  example  considered  below,  these  assumptions  are 
easily  verified  in  actual  applications. 

A  standard  technique  (see  e.g.,  [FI])  for  reformulating  (2.1)  as  an  abstract  Cauchy 
problem  is  essentially  to  set  A'  =  where  the  inner  product  on  A  satisfies 

(xi,x2)x  =  ^AI/2jq,  Ax/2x2^  and  to  take  the  state  space  to  be 

Z  =  X  x  V  x  L](-r.  0;.Y). 

The  approach  we  take  here  is  similar  but  it  allows  more  flexibility  in  the  choice  of 
state  space  in  that  X  is  not  necessarily  contained  in  Y,  but  X  will  be  in  one-to-one 
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correspondence  with  a  subspace  of  Y .  We  remark  that  explicit  knowledge  of  the 
square  root  of  A  is  not  required;  it  is  only  necessary  to  know  that  it  exists. 

With  this  discussion  in  mind,  let  5  be  a  subspace  of  Y  containing  P(.4),  and  let 
cr(-,-)  be  a  symmetric  bilinear  form  on  S  such  that  cr(yy,y2)  =  ^a.4y, ,  y2^  whenever 
y\  €  T)(A)  and  y2  6  5.  Let  A'  be  a  Hilbert  space  and  j  :  S  — >  A'  a  bijective  linear 
operator  such  that  j-1  :  A  — >  V  is  continuous  and  (x\,x2)x-  =  cr(j _1  T] , j_1.r2). 
Define  A  :  V(A)  C  X  —  Y  by 


V(A)  =  {j  6  X  |  j~'x  €  V{A)}  .  A  =  -aAj-K 


With  W  as  defined  above,  set  Z  =  X  x  Y  x  W.  Then  ;(/)  6  Z  satisfies 


where  A  is  given  by 


Jt  ~(0  =  -4-(0  +  col(0,/(0,0), 


V(A)  =  y  €  Z 


y  €  5,  w  e  V{D). 

(*  +  S-T9o{s)u'{s)ds^j  eT>{A) 


A  y  = 


j  47 

A  (*  +  f°rg0(s)u'(s)ds') 
jy  +  Dw 


THEOREM  2.5.  The  operator  A  generates  a  Co  semigroup  on  Z. 


PROOF:  Let  ,4,  D  and  j  be  as  above  and  consider  the  operators: 


V(C)  =  H7  Cw  =  J  g0(s)w(s)ds ; 


i  :  X  — >  U7  given  by  [?x](s)  =  x. 
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In  order  to  apply  Theorem  2.5,  we  must  show  that  A  can  be  factored  as  A  =  AqA). 


Since  A  :  V(A)  C  Y  — >  Y  is  positive  and  self-adjoint,  it  has  a  positive  scjuare  root 
A1/7 .  Define  .40  and  .4]  by 

V{A0)  =  V(A1'2),  A0  =  -aA'/2 
V{AX)  =  {. T  e  -V  I  j~lx  e  I >(^1/2)},  .4,  =  A 

Clearly,  .4  =  .40.4].  Set  C\  =  >4 1 C .  We  now  verify  the  conditions  of  Theorem  2.2. 

(1)  Clearly  i  is  continuous.  By  assumption,  j-1  :  7v(_/)  =  A’  — *  V  is  continuous. 

(2)  Since  5  2)  P(.4)  and  .4  is  densely  defined.  5  is  dense  in  V’.  Suppose  x  _L  D(A). 
Then  for  all  .r  E  P(,4). 


0  =  (x.x)  v  =  a(j  lx.j  '?)  =  a  (%'  , 

which  implies  that  ^.4y,;_1x^>  =  0  for  all  y  6  P(.4).  But  .4  is  self-adjoint  and 

one-to-one.  so  1Z(A)  is  dense  in  Y  ((Ft,  Theorem  13.11]).  Thus  =  0  which 

implies  that  x  =  0.  Therefore,  T>(A)  is  dense  in  A'.  Finally,  V[C)  fl  V(D)  = 

V(AC)  f~l  V(D)  =  //);(  — r.  0:  P(.4))  which  is  dense  in  IF. 

(3)  Obvious  since  7 Z(j)  =  X. 

(4)  We  already  know  that  D  is  closed.  To  show  that  F  is  closed,  let  j  €  V{F) 

for  n  =  1.2 . let  (*”  )  -*  and  f(^")  =  yn  -»  y  as  n  -»  oo.  Set 

i/n  =  ;_1  (in  +  I^r9a{s)Wn(s)ds^j  and  y  =  ;_1  (x  +  f°r  g0(s)u'(s)ds'j  .  Then 
yn  £  P(.4)  and  —aAyn  =  f(~  tn  )  ~  Vn  — >  y  as  7?  — ►  dc .  Since  .4  is  closed  and 


y  6  £>(-4)  and  —  ct Ay  =  y.  But  this  implies  that  x  +  gQ(s)iv(s)ds  6  ' D(A ) 
and  /I  ^x  +  f°r  ga(s)iu(s)ds j  =  y,  so  F  is  closed. 

(5)  Does  not  apply  since  7£((j2)  C  0. 


v 


U'{s),  n'(s))  ds  <  0 


.v 


bv  Lemma  2.4. 


(7)  Again  we  take  A0  =  1.  For  x  6  P(Lj),  (  j~Jx  1  e  T>(A),  and 

V(l-e6)x. 


rO 

L i.r  =  j_1x  —  .4  ^x  +  J  <7q(5)(1  —  es)xds  )  =  j_1x  —  qjA.t 
where  Qi  >  0  is  as  defined  above.  Define  T  by 


V(T)  =  {y  e  5  |  jy  e  P(Z,  )}.  T  =  I  +  QjQ-4. 

It  is  easy  to  see  that  T>(A)  =  P(T)  and  7Z(T)  =  'F(L\).  By  the  Cauchy-Schwarz 
Inequality,  for  y  6  V(T),  |jr?/|J  •  |jy|j  >  |  ( Ty,y )  |  =  ||y||2  +  QjQ  ^ Ay,y ^  >  ||y||2 
which  implies  || Jj  >  ||y||,  and  so  T  is  one-to-one.  Also.  Tm  =  T  so  TZ(T)  is 
dense  in  V. 

(8)  The  proof  that  (/  —  D)[D{C)  D  V(D )]  is  dense  in  IV  is  the  same  as  the  proof 
given  in  Section  2.1.  | 
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EXAMPLE.  A  viscoelastic  shaft  with  tip-mass. 

Consider  a  viscoelastic  shaft  of  length  /  fixed  at  one  end  and  with  a  tip-mass  at 
the  free  end.  The  equation  describing  the  motion  of  the  shaft  (see  [BMC]  where  the 

kernel  g(s)  is  assumed  to  be  in  //')  is 

q2  Q  '  Q  f  o  c) 

afa2yi't'X^  =  lh  Tfcy(t,x'i  +  J  +  s'x)ds  +K*)ux(th  (2.3) 

while  the  boundary  conditions  are  given  by 


?;(/.  0)  =  0.  (2,1) 

d2  r  d  f°  d  1 

An^i/(C/)  =  -  T—y(t.l)  +  J  g(s)  —  y(t  +  sj)ds  +  u2{t).  (2.5) 

Here  y  is  the  angular  displacement,  o  is  the  product  of  the  density  of  the  shaft  with 
its  polar  moment  of  inertia,  r  is  the  product  of  the  shear  modulus  and  the  polar 
moment  of  inertia.  Jm  is  the  moment  of  inertia  of  the  tip  mass,  and  the  delay  r  >  0 


is  assumed  to  be  finite.  Let  Y  =  R  x  L2(0J)  with 


Then  (2,3)  -  (2.5)  can  be  written  as 


y  +  A  Tv+  gi*)ysds  =  f(t) 


19 


Clearly,  T>(  A)  is  dense  in  Y.  Let  (  ^  "W  T>(A).  Then 


v(i)m 


/Vv  =  [  (v)'2 

Jo  Jo 


> 


21ml  +  °l2 


< 


>  0. 


Thus.  .4  is  positive  and  - - - - 

2  Jml  +  a  l2 

(and  .4'1  is  continuous).  Let  (  ^  )  ’  (  ^.2  )  €  72>(A).  Then 


so  .4  is  one-to-one 


.4  (  n 

fi 


'  r‘ 

>'l.  _  / 


■\{l)^2(l)  ~  /  L,  C2  =  /  c,v2 


so  .4  is  symmetric.  Let 


(l 

\v 


£  V.  Define  v(x)  — 


'  o  Ut 


av(0^s  +  1ml 


dt. 


Then  £  P(.4).  and  .4  J  =  J ,  so  71(A)  =  Y.  Therefore,  by 

Theorem  13.11  in  [R],  ,4‘  =  A.  Finally,  V(A~1)  =  71(A)  =  Y,  so  .4_1  is  closed. 


and  hence  .4  is  closed.  Now  let  5  = 

ri 


D « y 


xj  £  H1l( 0,/),  L-(/)  =  7f,  and 


define  a 


°  /  ^V'2  ^ 


(J,  )  ’  (  ?2  )  €  5-  Th(?n-  5  -  an<1 


<7(2/1- 1/2)  =  ( ctAyi.yij  whenever  r/j  €  2>(.4),  y2  £  5.  Let  X  =  H[(0,/)  with 
(xj.x2)a.  =  q  J^XjXj,  and  define  j  :  S  — >  .V  by  j  =  V’-  Clearly  j  is  a 

bijective  linear  operator,  and 


Ul-fdx  =  a  jo  •r'lX!  =  "  ((T/1  )  ’  (T^))  =  1,ri'J 

For  x  s  .V.  Ilr'xli;.  =  +  0  f'x1  <  (Ji^)  ||x||2v.  Tl.  -s.  J~'  is 


continuous. 
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Now,  define  Z  =  7/2(0,  /)  x  R  x  L‘( 0,/)  x  L](-r,  0;  7/2(0, /))  with 


The  operator  >4  generates  a  Co  semigroup  on  Z  by  Theorem  2.5. 

3.  Approximation.  In  this  section  we  consider  the  problem  of  finding  approximate 
solutions  to  equations  of  the  form  (1.1).  Let  S(t)  be  the  C0  semigroup  generated  by 
A.  We  construct  a  sequence  (Zn .  Pjf .  An)  where  Zn  is  a  finite  dimensional  subspace 
of  Z,  P%  is  the  orthogonal  projection  of  Z  onto  Zn,  and  An  generates  a  Co  semigroup 
Sn(t)  on  Zn.  We  then  show  that  Sn(t)P^z  — >  S[t)z  as  n  — ■*  oc  for  all  c  €  Z  using 
the  following  version  of  the  Trotter-Kato  theorem  which  follows  from  Theorem  4.2  in 
[P,  Chapter  3]. 

THEOREM  3.1.  Let  .4  e  G(M.rf)  be  the  infinitesimal  generator  of  a  C0  semigroup 
S(t)  on  a  Hilbert  space  Z.  For  n  —  1,2,...,  let  Zn  be  a  finite  dimensional  subspace 
of  Z  such  that  Pn  /  as  n  — ♦  oo  where  Pn  is  the  orthogonal  projection  of  Z  onto 
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ZJl.  Suppose 


(Hi)  .4'“  €  G '(M'S)  is  the  infinitesimal  generator  of  a  C0  semigroup  Sr,{t)  on  Zn  for 
n  =  1.2 . and 

(H2)  for  all  z  £  Z.  (XI  —  An)~iPnz  — >  (XI  —  -4)_1c  as  72  — >  oc. 

Then  for  all  z  £  Z.  Sn(t)Pnz  —>  S(t)z  as  n  —>  oc,  and  the  convergence  is  uniform  on 
bounded  t-intervals. 

We  would  like  to  construct  a  convergent  approximation  scheme  for  the  thermo¬ 
viscoelastic  system  considered  in  Section  2.2.  Nevertheless,  for  the  approximation 
scheme  and  convergence  proof  presented  below,  we  restrict  our  attention  to  the  ab¬ 
stract  viscoelastic  system  (equation  (1.1))  for  which  the  operator  A  is  given  by  (2.2). 
The  convergence  proof  for  the  complete  thermo-viscoelastic  system  involves  a  modifi¬ 
cation  of  the  proof  we  give  here  and  can  be  found  in  the  thesis  [L].  This  modification  is 
rather  technical,  and  yet  is  a  straight  forward  extension  of  the  proof  we  present  below. 
Therefore,  in  order  to  conserve  space  we  present  only  the  proof  for  the  viscoelastic 
model.  If  we  define  the  operator  -40  :  T>(A0)  C  X  x  Y  —*  X  x  Y  by 

I’M,,)  =  j(J)  €  -V  x  y  |  X  €  Z>(.4),  yesj, 

then  we  can  write  A  in  the  form 
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which  suggests  a  two-stage  approximation  of  A'.  We  first  approximate  .4o  bv  discretiz¬ 
ing  the  spatial  variable,  typically  by  means  of  spline  functions.  We  then  approximate 
—  bv  discretizing  the  delav  variable.  In  this  paper  we  will  use  an  averaging  scheme 

OS  - 

for  the  second  stage.  We  follow  the  construction  given  in  [Flj  except  that  we  do  not 
require  a  uniform  partitioning  of  the  interval  [  — r,  0]. 

Let  us  now  proceed  with  the  first  stage  of  the  approximation.  Define  the  bilinear 


form  <70(-.  • )  on  X  x  5  by 


<70 


where  a  is  the  bilinear  form  on  5  discussed  in  Section  2.3.  Observe  that  for  j  £ 
D(.40)  and  (H)  e  X  x  5, 

Now  for  each  positive  integer  N.  let  Xs  and  V  'v  be  finite  dimensional  subspaces  of  .V 
and  Y  with  V’A  C  5,  and  define  IVA  =  L2g(  —  r,  0;  A  A  ).  We  define  :  A  A  x  V  A  —► 
A’ A  x  by  restricting  <7q  to  .Y  A  x  V  A ;  i.e., 


«/./) A.xV-  =  (to  (/./)  for  u‘v, t,A'  €  Xs  x  Ys 


Now  set  Z‘v  =  A  A  x  V"  v  x  U’A  and  define  A A  :  P(^4A)  C  ZA  — ♦  ZA  by 


V(AS)  =  XN  x  Yn  x  HlR{-r.  0:  Aa'). 


-4*  (  yv 


( 


4V 


+  /  (Ja{s)  ^  U'  )  ds 

■  \  ,  d  ,v 
J»  +^U' 


\ 

) 


23 


For  each  positive  integer  M  partition  the  interval  [— r,  0]  into  subintervals  [t1^ ,  J 


j  =  1,2, ... .  M ,  where 

— r  =  t\\  <  <  •  ■  •  <  t\*  <  Iq  —  0.  (3.1) 

We  will  sav  more  later  about  how  the  /A/  are  chosen.  Set  a  A^  —  ,  —  t*1  for  j  = 

1,2....,.!/,  let  y  A/  denote  the  characteristic  function  of  for  j  =  2 ...... U, 

and  let  \A/  denote  the  characteristic  function  of  [fA/,0].  Let  D^(t).  i  =  0.1, . M 

be  the  usual  linear  spline  functions  satisfying  Z?A/(/A/)  =  6,j.  Define  the  finite  dimen¬ 
sional  subspaces  H’AA/  and  1FV'A/  of  IF  by 

WS'M  =  |u:  €  w  |  u-  =  £  6  -VA  }  , 

If  v,w  =  jut  €  IF  iv  =  fj  6 ;VfiA/,  e  A'  v  j  . 

Define  5A’A/  :  -  WNM  by  D*""wNM  =  E(=i  ptr  (4&  -  if)  \f  where 

A/ 

rt',V’A/  =  6Af  j9a/  and  b f  =  0.  Define  the  isomorphism  fA'A/  :  l'FA,Af  — +  IFA,A* 

i=i 

by  =  Ziii  bi*Xi*-  Now  define  DA'Af  :  1FA'-A/  —  1-F‘V'A/  by  £>A’-W  = 

DA"U  (?A  A/)  \  To  complete  the  approximation,  set  ZA,A/  =  .VA  x  V'A  x  H/A,u. 
and  for  rA,A/  =  (xA  ,  yA  ,  u,A,Af  )T  £  ZA,A/,  define 
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Af  =  {j6AM|af  <r,W-‘,+c>J/2}. 

If  j  £  AA/,  then  ( ga )A/  <  — ,  and  contains  at  most  3/1_£j  elements  of  AM. 

Furthermore,  oA^1(^70)A,  <  (ga) Al,aA/  for  J  —  2,3 . U.  and  if  j  €  Aff,  then 

1,2 . j  -  1  €  A”. 

EXAMPLE  3.2.  Suppose  t*1  =  for  J  —  0. 1 . M.  Then  o;'f  =  —  for  all 

JM  1  M 
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j  6  A",  so  (A2)  is  satisfied  with  (\  =  1,  (2,C  >  0  arbitrary  since  A**  =  0  and 

Example  3.3.  Let  C  =  f^r  ga(s)ds,  and  suppose  that  the  l™  are  chosen  so  that 

c 

(go)*1  =  —  for  j  =  1.2 . M.  Then  it  is  easy  to  check  that  for  fi  =  <2  =  5  the 

partition  0 A/  satisfies  (A2)  for  all  positive  integers  M . 


We  will  refer  to  these  partitions  as  the  “uniform  mesh”  and  the  "non-uniform  mesh,'' 
respectively.  The  non-uniform  mesh  of  Example  3.3  is  the  partition  suggested  by 
Fabiano  and  Ito  in  [FI].  The  convergence  proof  we  give  below  is  a  modification  of  the 
proof  given  in  [FI].  The  major  difference  is  that  we  must  handle  each  estimate  in  two 
parts:  one  where  the  length  of  the  interval  is  small  (in  this  case  our  argument  is  the 
same  as  that  of  Fabiano  and  Ito),  and  the  other  where  the  integral  of  g  is  small. 

Let  P^'*1  denote  the  orthogonal  projection  of  Z  onto  ZSS1 .  In  order  to  apply 
Theorem  3.1.  we  must  show  that  Pz  "  *  lz  as  .V,  M  — >  oc.  For  r  =  ( x.y .  w)1  €  Z , 

Pz’*1  z  —  [p^  x.  Py  y.  P^'*1  trj  .  where  P^'*1  is  the  orthogonal  projection  of  IF  onto 
IFAA/.  Since  we  assume  (Al)  it  is  sufficient  to  show  that  P^'*1  I\y. 


LEMMA  3.4.  For  all  h  6  H  .  P^.'*1  h  — »  h  as  .X.M  — *  oc. 

PROOF:  For  w  6  IF  and  t  £  [ — r.  0] .  set  u;A(f)  =  P\tw(t).  and  define  w1*'*1  by 

vSM(t)  =  E^V(L-,)\iU(0.  For  ic  e  V(D2) 

1=1 


,y„\/  ^  ,v  .  .v  X.M 

iv  -  ir  1(,  <  iv  —  w  u.  +  ir  -iv  |lv. 


where  the  first  term  on  the  right-hand  side  tends  to  zero  by  the  Dominated  Con- 
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vergence  Theorem.  We  write  the  second  term  as  || u,A  —  u,,VAf|j^.  =  S\  +  S?  where 
S-,  —  ||?cA  (s)  —  w's  )  j|  'v  (Is  for  j  =  1.2.  Following  Fabiano  and 


Ito  (but  replacing  r/M.  the  length  of  each  interval  using  the  uniform  mesh,  by 

r/.\/F  +  £l)/‘)  we  gC( 


S\  <  r\  D 


2  * 

w  n. 


’  3  (.\/F  +  t>)/‘) 


0  as  M  —*  dc. 


Now.  for  s  £  [— 7\ 0]. 


I 5) || _v  <  ||u-(5)||i 


/  =  /  “7=777 

b  .V  ||  A  \/<7o(0  |_v 

-d=4e  fs.Ki  l|27u-(Olli  <7f  <  -f-  IIDu-lli. 

V<7a(sc)  Js  9a(-r) 


Thus. 


S,  <  2 


&>($)  II-a(^)|Ia  +  II-a(^1)IIa  * 


4r||Du’||(r  ^  f‘-'  ■ irC 

—j - 7~  2_.  /  9o(s)ds  <  — — 


0  as  3/ 


Hence,  for  te  G  V[D2).  iv  -  P^'AI  ic  <  ||tr  —  u,A'A/||u.  +  (io  —  reN''U)  < 

2  Iju'  —  u,*A -a/ 1 1  — >  0  as  .V.  M  — ♦  cc.  Since  V{D 2)  is  dense  in  IF  and  II  Jj  <  1. 


h  -  P™h 


0  as  .\.M  —>  dc  for  all  h  G  IF. 


From  Theorem  2.5  we  know  that  A  G  G'(l.O).  In  order  to  show  that  >FN ,M  g  (7(1,0) 


for  all  N,M,  it  is  sufficient  to  show  that  _4n,a/  js  <J issj pat i \e  in  since  ZN"U  is 
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(  X  \  M 

finite  dimensional.  Let  zA,A/  =  I  yA  with  u'A,A*  =  u>Af \A/ .  Then 

V  /  -=i 


\N,M  „N,M  _.V.Af 


)z  =  a° 


*"  +  £(*>)«'  fxA 


’  h/N 


A/  M  \]  A/ 

in  -  w:  \  - 


+  ».w  +  »?v 


A/ 

E1  ,  \M  /  At  At  M  \ 

~Jl(9o),  “  U’t-  ,u-t  )  v 


A/  ,  ,jv/  r 

^  V-'  \9a  )i  II  At  II  II  At  II  II  At 1 1 2 

-  /  .  At  Ila>-1  II.Y  ‘  llU'«'  II.Y  IK  II  A’ 

li  '  t- 


4z^[M-ihuiiL 


i  r^-1 

=  ;  likx 

L  i=i 


At  ||2  f  (da)  i  +  i  (<?cr)f' 


At  2  \9a)\f  ^ 

.y-ta r  s 


where  we  used  the  Cauchy-Schwarz  Inequality  and  the  inequality  2 ab  <  a2  +  b2.  and 


from  (A2)  the  fact  that  (^o)f+i/Qf+i  <  (</o)fVQi'/  f°r  ?  =  1,2,...,  3/  —  1.  Thus,  we 


have  established  the  following  result. 


Lemma  3.5.  If  A**1  is  as  defined  above,  then  .4A,A/  £  G(I,U]  lui  all  .V,  \f . 


All  that  remains  to  be  done  in  order  to  establish  convergence  is  to  show  that  (A/ 


^iV,A/)-ipA’.A/,  for  a|]  ;  £  Z.  Lor  ReA  >  0  and  c  =  (  i/>  ,  consider 
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the  equation  {XI  -  A)' 


y  I ,  or  equivalently, 

U’  / 


Ax  -  jy  =  p, 
Xy  —  A  ( x  -h 


J  9o(s)w 


{s)ds  = 


\W-Jy-CXL  =  I,. 


From  (3.4).  u'{s)  —  tx^~^{jy  +  Ms ))d£,  and  from  (3.2).  jy  =  Ax  —  r*.  or  y  = 
Aj-1x  —  V-  Substituting  into  (3.3)  and  using  a \  as  defined  above  and  the  fact  that 
Xexls~^d£  =  1  —  eA*  we  obtain 

A2j-1x  -  .4  ^q,\x  -  J  ga($){XI  -  D)~}  [y  -  h{s)]ds^j  =  v  +  A j'V-  (3.5) 


If  we  define  A(A)  =  A 2j  1  —  —  A  ^ r  +  f°r  exsg(s)dsj  =  X 2j  1  —  a.\.4.  then  for  x  6 
V(A), 


A(A)x  =  y  +  A  j 


i  'y  —  A  J  gQ 


[s){XI  -  D)~l[p  -  h(s)}ds. 


Now,  let  P7 


V'V  J ,  and  consider  the  equation  (A/  —  <4A,A/)  P 
hXM  ) 


.wuy1  pA'.At,  _ 


I  .v  /  X 


Vs’ 

by  (,4'N  x  N  .  y' 

'  >v  =  U 

jys 
4  vx 

■OforU;'b 

)  6  A  *'  x  )■ 

the  equations 


\zx-jyN  =  'ps, 


M 


Xy* -A*  l  x"  +  £>)«'  =v*A\ 


wj 


Au-f 

.  V  1 

-,y 

(™i-i 

M 

and  h' 

M 

=  E«* 

\\M  _ 

EW 

1=1 

i=i 

(a+7? 

) 

=  —  mA/ 
q;VU- 

M. 


(3.C) 

(3.7) 

(3.8) 


or 


■|'U  -  (1  +  q;UA)  [u*;Yj  +  a,u  (jyA  +  /i;u)]  for  ?  =  1.2 . U. 


,-i! 


17  0  +  q/VA)  J  aA/  (ji/A  +  From 


where  u'0u  =  0.  By  induction.  irA/  =  £ 

k=l 

(3.6).  jyA  =  AxA  -  ^A  which  implies  yv  =  \j~lxN  -  j~'^s .  Substituting  into  (3.7) 
we  obtain 


M 


#r'*S-A*  U'V  +  A^:(,0);^ 


i=i 

M 


k=  1  U=* 


n  (1  +  af'A) 


Q<7X  V 


p-v  +  A;-v-v-^(^);u^vE 


1=1 


A=1  L/=t 


IJ  (1  +  a;v  A) 


Q 


ns-h").  (3.9) 


By  induction,  A  ^ 


I7(l+af'A) 


M  W~l 


k=  1  L/=* 

Aa’-u(A)  =  A2;-1  -  ,4V 


Q‘k  —  1  —  17  0  +  q7* A)  '.  Define 


k-\ 


M 


1 + ^(soYi1  1  -  n  o + q^a) 


i=i 


A  =  1 


Then,  since  j2(g0)-’  =  --  [  g(s)ds  =  -  -  1. 
.=  i  a  J- T  a 


Aaaw(A)  =  A2;-’  -  --AA 
a 


T  + 


j  g{s)e-''{\.s)ds 
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(3.10) 


where 


^%#(AlS)  =  f;(n(l+o"A)-,)x;V 

1=1  ' k~i  ' 

Let  (XI  -  Ds'M)~l  (v?'v  -  hN-")  =  jr  &MV;V.  Then 


M 


M  r 


\c M _ L_  (c M  cM 

a" 


(e,-^)  \.A/ 


1=1  1=1 

which  implies  ;pA  —  /)A/  =  A{  A/  —  —  £A/)  for  2  =  1.2 . A/.  Thus.  <fA/  = 

i 

(l  +  AaA/)  1  [s^li  +  Q,W  (-pA  —  /fA/)],  or.  by  induction. 


CM 

s. 


=  E 

k  —  l 


ri  (i + xqi‘) 


:=k 


A/  /  „.V  i  A/\ 

Qk  (r  "k  ) ■ 


Now.  f  ga(s)  E^;U\;V  =  E(</o)i'/^,A/,  so 

2 - r  1=1  1=1 

j“  sM  (M  -  Ds  >’)~'  * 


A/ 


E^E 


1=1 


A=1 


A/\-> 


n  (i + wo 


,/=t 


a?  (*"  -  h?) .  (3.11) 


Therefore,  from  (3.9).  (3.10)  and  (3.11)  we  obtain 

Av,a,S  A  )xv  =  L’V  +  A  rVV  -  AX  J\o(s)  (A/  -  (v*‘V  -  ^'V-M)  *•  (3-12) 

In  order  to  complete  the  proof  of  convergence,  we  must  show  that  xjV  — ►  x,  yA  — >  y 
and  u>'v,A/  — >  re  as  .V.  M  — *  oc.  First  we  need  the  following  lemmas. 

Lemma  3.6.  For  A  >  0.  (XI  -  T>V  A/)_1  PAA'/i  -  (XI  -  D)~'h  for  all  h  €  IF. 

PROOF:  Let  w  =  (XI  -  D)-]h  and  wXM  =  (XI  -  DNM)~l  P„M h.  Set  u-‘VlM  = 
(P'v,a/)-,u,a.a/.  ancj  ^ a he  u’A,A/  =  P,AA/ (Dw)d£.  It  now  follows  that 


i™  («)-v-A/  -  CrVA/)||lv  <  A 
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P^M  (u-  - 


IV 


(3.13) 


Now, 


IV  —  1‘  '"it'  ”  ...  <  u>  —  U'  '  ...  +  u>  —  ?  u-  u. 


w  lie  re 


-..V,A/!|2 


Die  -  Pu""(Dw)\  J‘ 


<  ~  Dw  -  Pw’m(Dw) j  *  —  0  as  .V.  M  -»  oc 


.  w 

hv  Lemma  3.1.  If  we  set  (Dir)  =  ^ien  we  Set 


||u-v-A/  -  rv"wu- 


-  E  /„  *wrt  (*  -  ',’T*  ■ 


.S  j  -f"  2  • 


To  estimate  the  term  5;  we  again  follow  Fabiano  and  Ito  (replacing  the  term  r/M 

o 

by  r/.U(I  +  f,,/2)  to  get  5,  <  u'1+-  ||Die||2r  -»  0  as  M  -»  oc.  Define  the  norm  |H|j 

on  IF  by  || i4*|| i  =  £|k(,)||i*.  Then  ||/tf"''(Die)j|*  =  ^  ||Ei=i  *  = 

m  /■*:-» . o  D.  . 2 


£  /"’  lk"ir,*  =  Eoi“lk.”l!V  f'°' 

.=  1  JfM  .=  1 

I 

imi?  =  J  ikwiii  (/-s  ^  ~ j  J  9M  i 


|ir(s)l|2v  c/s  =  — t - r  ||u'||(r 

Fo(-r) 


and  so.  for  all  j , 


>f  IldMIf  <  ||^:j'(d«-)||'  <  -F-  ||D-’'(o-)||l  <  ^F_  nDu'iiJ,. 


'  9o(~r) 


In-  ga[-r) 


Thus,  for  5  €  [1**  .t*h),  ||£;U  ||.v  (-s  ~  lj' )  -  (Qf)  ll^ll.v  -  ^_r) 


Hence, 


v-  r  ,  7'  II Du’ ||  m-  C  n 

S  E  ^  /„  *<*>*  £  Jpitk  - 0  “  " 


and  so  \\w  —  I 


0  as  A’,  M  — *  oo.  Thus,  from  (3.13)  it  follows  that 


trVM  —  i N — *  0  as  A\  .U  — *  oc.  Therefore, 


V,A/  ^  .V.  W  -  A’,  A/ 

rr  —  u’  lv  <  ir  —  2  1 


+  |rVA/u-VA/  —  U"VM  lv  — >  0  as  .V,  A/  — *  oc. 


Lemma  3.7.  For  A  >  0.  J  gQ{s)  |eAS  -  eA,(  A.  s)|  ds  —  0  as-  .U  —  oc. 


PROOF:  By  definition  of  ru(A.>). 


<?c (-^ )  eA"s  -  ts1[ A 


.  A  •  -<  )  |  d.s  =  ^  e  A"  —  JT  (  1  +  QjW'M  ^ 

tr  A"  f=i 


=  5i  +  S2.  Let  c  >  0,  and  choose  A/0  large  enough  that  if  M  >  M0 ,  then  y—  <  t, 

-  l  <  1  and  <  ‘  For  >  ,Uo. 

2  2 .  w c  1  2 


Si  <2^2  [  9°(*)d 


-  M e* 


Let  Aa/  =  {1.2, ...  .n)  and  suppose  s  €  [fA/,0].  Then  for  some  ?  6  AA/,  txl  <  s  < 
which  implies  that  eXt™  <  exs  <  eA(*-i.  so  0  <  exs  —  eA(<Af  <  eXt*'  ^eAo>W  —  1^  <  ^. 
But.  txl  =  —  ]T^_,  Q j1  which  implies 

.‘•.--no+nfA)-  =  ne-'-“  -  n  (i + m'^)" 

j=i  j=i  r=i 

<  n(i+o*'A)-ne-r , 


Now.  t  =  (1  +  Ao  J1)  +  -  f  { A  o  * '  \ 1  for  some  f  ,  between  0  and  a**.  Set  a }  = 
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(l  +  Aaf)  and  F,  =  -cXi]  (A nA/)2.  Then 


ri(i+0?'A)-n^  = 


b  =  i  j=i 


By  induction,  FK  -  n  K  +  M  =  -  E  bj  n  (<U-  +  M  n  a‘-  Thus, 

J  =  l  J  =  1  jsl  A-J  +  l  /=1 


noTo^-n^'  =  (Aa^r  n  eXa* 


Jt=J+l  /=) 


£F',Z(An”)T^'r-  M 


Xr  \2  eAr(  Ar)2  ^  c 
\f{ i  +  c,)/2  )  ~~  2.\/e'  <  2 


since  e-^  fl  f'a‘  fiV^A)  <e'  ''  fl  e  W'  fl  eAo"  =  eA‘"  <  eAr.  Hence. 


k=}+ 1  /=1 


*=J+1  /=1 


■?i  <  X  j'~  9*(s)  ^  -  C'W;M  +  eA<i'f  -  n  (]  +  q;/A)_1  * 

■  sa*'  ■  L  j=i 

A/  /-o 


21  /-O 

<  «  Y]  /  ga{s)ds  =  t  /  #>($)*. 

trAw 


We  are  now  read)’  to  prove  the  main  result  of  this  section. 


THEOREM  3.S.  For  all  :  £  Z,  c4‘  M(.P?N"',~  — 1 1  S(t):  as  .\,M  — >  oc.  uniformly 
bounded  t-intervals. 


PROOF:  As  remarked  above,  we  have  only  to  establish  ( H 2 )  of  Theorem  3.1.  Fur¬ 
thermore,  it  is  sufficient  to  show  that  x,N  — *  x.  If  we  define  the  bilinear  forms 


//(•.■)  on  A  by  /i(X].j2)  =  A 2  (j  1  x  j .  j  ’x2)r  +  a.\  (x  i-x2)v  (it  is  easy  to  verify 


that  ft(x i.x2)  =  (_A(  A )j, .  x2)y  if  X]  6  P(.-l)).  and  //A/(-.-)  on  A  N  by  /rA,(x i.x2)  = 


(A‘v,a/(A )xuj  1 0'2)1,  ,  then  for  xi,x2  G  -VA, 

|/ia,(xi,x2)  -  /i(x,,x2)|  <  ^ J  ga{s) \exs  -  eA/(A,s)|  ds'j  ||xi||  v  •  ||x2||  v  .  (3.14) 

and  for  x  6  Xs , 

/iA,(r,  j)  =  A2  ||j~1x||y  +  a  a  H-rllx  ■  (3-15) 

By  (3.5),  we  have  for  all  u  G  -V, 


35 


from  (3.15)  the  fact  that  a\  xA  —  xN  y  <  //M(5'N  —  xN,zs  —  xN),  we  get 


i.v  ,..vi 


I.Y 


< 


a  .\ 


(^2  I|i  'll  +  oA/^  ||xA  —  x||  v  +  9ois)  \(  SS  -  eM(^"s)|^  11*11 


+  \\j  1 1|  •  11^  -  t,A  ||v- +  a  ||j  ‘|[  -||* -</!!*  +  (/  9°(s 


)ds 


1/2 


(A/  -  -  /,)  -  (A/  -  Dy',rlP^,(^  -  h) 


irJ 


as  .V.  M  — ►  oc  by  Lemmas  3.(i  and  3.7.  so  J|.r  —  ,r  N  ||  v  — »  0  as  A .  M  — +  oc.  | 

We  are  interested  in  applying  this  approximation  scheme  to  the  optimal  LQR  prob¬ 
lem.  but  it  is  well  known  that  convergence  of  the  forward  problem  ( i.e. ,  Theorem 
3.8)  is  not  sufficient  to  ensure  convergence  of  the  gain  functionals  (see  [G],  [BIP]). 
Convergence  of  the  adjoint  semigroups  as  well  as  other  properties  (uniform  stabiliz- 
abilitv.  etc.)  play  a  central  role  in  the  development  of  convergent  methods  for  LQR 
problems.  If  g(-)  €  L2{—  r.O)  then  Z  can  be  re-normed  by  using  the  weight  e(s)  =  1 
in  place  of  g(s).  If  Ze  denotes  the  resulting  equivalent  space,  one  can  establish  the 
following  adjoint  convergence  (see  [M]  for  a  proof). 

Theorem  3.9.  If  g(s)  G  Z2(-r.  0).  then  for  each  :  G  Zc. 


eK*T  'py"- 


S“(t)z  as  A’,  M  — ►  co. 


Moreover,  this  convergence  is  uniform  on  hounded  /-intervals. 

Recently.  K.  Ito  has  announced  [I]  a  proof  of  this  result  that  does  not  require  the 
additional  assumption  that  <j  ( ■ )  G  L2(  —  r,  0).  In  general,  the  question  of  preserration 
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of  stabilizability  uniformly  under  these  approximations  is  not  answered.  However,  for 
certain  special  forms  of  g(-)  and  for  problems  with  additional  damping  terms  (Kelvin- 
Voigt),  one  can  establish  uniform  stabilizability.  More  will  be  said  about  this  problem 
in  the  next  section. 


4.  Numerical  Results.  We  turn  now  to  an  optimal  control  problem  governed  by 
the  basic  thermo-elastic  equations  with  viscoelastic  and  Kelvin- Voigt  damping  terms. 


For  a  rod  of  length  1  the  equations  of  motion  become 

d2  ,  d 

r’Wy(Ul)  =  TI 


(A  +  2,0  +  i  j +  s.xHs  - 


Q 

-  q(3A  +  2 fi)-O(t.x)  +  b{x)u{t)  (4.1) 

ox 

d  c)~  d* 

pc—O^-x)  =  K^pO(t,x)  -  ^oQ(3A  4-  2g)-^-y[t,x)  (4.2) 

. 

where  for  0  <  3  <  +ec.  the  teim  3  y(t.x)  provides  Kelvin- Voigt  damping.  If 

otox1 

3  =  0  and  g(s)  =  0.  then  (4.1)  -  (4.2)  become  the  classical  equations  of  thermo- 
elasticity.  If  g(s)  7^  0  satisifies  Hypothesis  2.3,  then  the  integral  term  provides  a  type 
of  “viscoelastic  damping'-  to  the  system. 

Gibson.  Rosen  and  Tau  [GRT]  considered  the  problem  with  g(s)  =  0  and  with 


Dirichlet  boundary  conditions  on  displacement  and  Neumann  boundary  conditions 
on  the  temperature.  In  this  case  it  can  be  shown  that  the  thermo-elastic  model  (i.e.. 
3  =  0  and  g{s)  =  0)  has  zero  as  an  exponentially  stable  equilibrium,  provided  of 
course  that  one  subtracts  out  the  constant  temperature  distribution.  Moreover,  the 
elastic  and  thermal  modes  decouple,  and  one  can  use  modal  expansions  in  analyzing 


37 


this  system.  Therefore,  if  3  =  0  and  r/(s)  =  0,  then  (4.1)  -  (4.2)  with  Dirichlet- 


Neumann  boundary  conditions  is  exponentially  stabilizable  (again,  after  "subtract¬ 
ing"  the  constant  temperature  states),  and  the  standard  LQR  control  problem  has  a 
unique  solution.  The  same  clearly  holds  for  the  corresponding  systems  with  ^  ^  0 
and  g(s)  ^  0. 

We  shall  consider  the  more  complex  problem  governed  by  (4.1 )  (4.2)  with  Dirichlet 

boundary  conditions  on  displacement  and  temperature.  In  particular,  we  impose  the 
boundary  conditions 

y(t-0)  =  g((.  1)  =  0  =  0(t,0)  =  0(t  A).  (4.3) 

If  3  >  0,  then  (4.1)  -  (4.3)  is  exponentially  stable.  If  3  =  0  and  g(s)  satisfies  the 
basic  assumptions  in  Section  2  with. 

1 

c0  =  (A  +  2/i)  +  -J  g{s)ds  >  0,  (4.4) 

then  it  is  known  that  (4.1)  -  (4.3)  has  zero  as  a  globally  asymptotically  stable  equi¬ 
librium  (see  [W.pp.  203-210]).  It  is  important  to  observe  that,  in  general,  results  of 
this  type  do  not  extend  to  problems  in  two  or  more  space  dimensions.  If  3  =  0  then 
it  is  still  not  known  if  (4.1)  -  (4.3)  is  exponentially  stable.  There  are  positive  results 
for  infinite  delay  problems  with  completely  monotone  kernels  g(-)  (see  [HW]). 

Recently  Hanson  [H]  has  shown  that  for  ,3  =  0  and  g(s)  —  0,  the  open  loop 
eigenvalues  of  (4.1)  -  (4.2)  are  bounded  away  from  the  imaginary  axis.  However,  the 
eigenfunctions  for  this  problem  do  not  form  a  Riesz  basis  for  the  natural  state  space. 


and  hence  the  question  of  exponential  stability  remains  unresolved. 

These  theoretical  issues  can  have  considerable  impact  on  computational  algorithms 
and  need  to  be  addressed  in  order  to  complete  the  convergence  theory  for  control 
applications.  Some  progress  has  been  made  for  the  thermo- viscoelastic  problem  with 
kernels  of  the  form 

t/f-")  =  -tqs/py^.  (4.5) 

for  q  >  0.  p  >  0  (see  [Lji.  although  no  results  exist  for  the  classical  thermo-elastic 
model  3  —  0.  <y(s)  =  0.  In  summary,  if  3  >  0  or  </(•)  has  the  form  (4.5)  with 
(4.4)  satisfied,  then  (4.1)  -  (4.3)  is  exponentially  stable  (hence  stabilizable).  and  if 
3  —  0  =  g{$).  it  is  not  known  if  (4.1)  -  (1.3)  is  stabilizable. 

We  shall  present  several  numerical  experiments  for  LQR  control  of  (4.1)  -  (4.3) 
with  various  values  of  3  and  b  when  y(-)  is  given  by  (4.5).  The  numerical  constants 
are  chosen  to  be  the  same  as  used  in  GRTj  for  an  aluminum  rod  of  length  1.  In 
particular. 

p  —  9.S2  x  10-2.  A  =  2.06  x  10"1. 

p  =  1.11  x  10"1.  q  =  1.29  x  10“3. 

c  =  5.40xl0-1.  k  =  7.02  x  1CT' . 

0o  =  GS.  q  =  30. 

p  =  1 .  r  =  1 . 

Observe  that  for  these  values,  the  constant  y  =  o(3A  +  2//)  takes  the  value  l.OSSx  10~3, 
and  lienee  the  coupling  between  (2.1)  -  (2.2)  is  "small."  Note  also  that  the  kernel 


If  g(s)  =  0  and  3  >  0,  then  again  there  exists  a  unique  optimal  control  law  in 
feedback  form  (4.9)  with  g(s)  =  0  =  J\4(s,x).  If  j  =  0  and  g(s)  =  0.  then  the 
most  one  can  say  is  that  if  an  optimal  control  exists,  then  it  will  have  the  form  (4.9) 
with  g(s)  =  0  =  I\A(s.x).  The  existence  of  the  optimal  control  has  not  yet  been 
established.  We  had  hoped  that  the  numerical  results  presented  below  would  shed 
some  light  on  this  question.  However,  as  will  become  evident  from  these  numerical 
experiments,  we  seem  to  have  raised  more  questions  than  we  have  answered. 

The  objective  of  the  computational  scheme  is  to  produce  numerical  approximations 
of  the  optimal  functional  gains  A'j(-).  A’2(-)>  A'3(-)  and  hi  particular,  the 

general  idea  is  to  introduce  an  approximation  method  (such  as  in  Section  3).  and 
then  to  use  this  method  to  compute  approximating  (sub-optimal)  gains  A'jN  (•),  K*  (■). 
A'jN(-)  and  A'^N  (-,-),  The  basic  questions  are  i)  Do  the  optimal  gains  exist?  and  ?7) 
Does  the  particular  approximation  scheme  lead  to  convergence  of  the  sub-optimal 
gains;  i.e.,  does  A’,N ( • )  — ►  A'.(-).  i  =  1.2, 3, 4  as  .Y  — ►  oc? 

The  computational  results  presented  below  are  based  on  the  numerical  approxima¬ 
tion  scheme  presented  in  Section  3.  We  used  the  non-uniform  mesh  as  described  by 
Example  3.3  for  the  thermo-viscoelastic  model.  In  particular,  all  runs  presented  below 
were  based  on  M  =  S  subdivisions  of  [—1,0],  and  the  corresponding  approximating 
finite  dimensional  system  becomes 


=y{t)  =  As*z*(t)  +  Bxsu(t) 


(4.10) 
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where  =  P^'SB  and  B  :  R1  — >  Z  is  defined  by  Bu  —  (0,  Ii(t)ii,  0,  0)7".  The  results 
are  the  same  for  M  =  1G  sub-divisions  of  [—1,0]  (a  nice  feature  of  the  non-uniform 
mesh  algorithm). 

In  the  cases  where  i)  holds,  Theorems  3.S  and  3.9  imply  convergence  of  the  func¬ 
tional  gains  provided  one  can  establish  that  this  scheme  preserves  stabilizability  uni¬ 
formly  in  A’  (see  [G]).  For  the  thermo- viscoelastic  model  with  </(•)  given  by  (4.5)  and 
the  thermo-elastic  model  with  3  >  0.  preservation  of  stabilizability  can  be  shown. 

In  order  to  see  the  effect  of  the  damping  models  on  the  open-loop  system,  in  Figure 
1  we  plot  the  open-loop  poles  for  the  first  8  clastic  modes  using  only  thermal  damping 
(i.e.,  3  —  0.  g(s)  =  0),  thermo-viscoelastic  damping  (3  =  0,  g(. s)  given  by  (4.5)  and 
b  =  6.0  x  104)  and  Kelvin- Voigt  damping  (3  =  l.S  x  10~'),  respectively.  Observe 
that  the  damping  is  extremely  small  in  all  three  cases.  On  the  other  hand.  Table  1 
shows  that  at  the  higher  frequencies  the  damping  provided  by  the  thermo-viscoelastic 
model  is  an  order  of  magnitude  more  than  the  thermo-elastic  model,  and  the  damping 
predicted  by  the  Kelvin-Voigt  model  is  two  orders  of  magnitude  greater  than  the 
thermo-elastic  model. 
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elastic 

Mode 

Thermo- Elastic 

3  =  0  =  g(s) 

Thermo- Yiscoelast  ic 
b  =  6.0  x  104 

Kelvin- Voigt 

3  =  l.S  x  10"7 

1 

—0.0000054  +  6.5S? 

-0.0000091  +  6.58/ 

-0.0000093  +  6.58/ 

2 

-0.0000060  +  13.17? 

-0.00001S4  +  13.17/ 

-0.0000216  +  13.17/ 

3 

— 0 .000007 1  +  19.79/ 

-0.0000325  +  19.79/ 

-0.0000423  +  19.79/ 

-1 

— 0.00000S6  +  26.-17? 

-0.0000501  +  26.47/ 

-0.0000715  +  26.47/ 

5 

— 0.00001 06  4  33.20? 

-0.0000699  +  33.20/ 

-0.0001095  +  33.20/ 

6 

— 0.0000129  f  40.02/ 

-0.0000905  +  40.02/ 

-0.0001566  +  40.02/ 

7 

-0.0000157  -f  46.93/ 

-0.0001113  +  46.93/ 

-0.0002133  +  46.93/ 

8 

-0.0000189  +  53.95/ 

-0.0001317  +  53.95/ 

-0.0002800  +  53.95/ 

9 

— 0.000022-4  +  61.10/ 

-0.0001514  +  61.10/ 

-0.0003574  +  61.10/ 

10 

-0.0000263  +  68.39/ 

-0.0001705  +  68.39/ 

-0.0004460  +  6  S .  3  9  ? 

15 

-0.0000491  +  107.5/ 

—0.0002565  +  1 07.5/ 

— 0.00 1 0S66  +  107.5/ 

20 

-0.0000662  +  151.5/ 

— 0.0003257  +  151.5/ 

-0.0021294  +  151.5/ 

TABLE  1.  OPEN  LOOP  POLES;  SMALL  DAMPING 


We  also  considered  cases  where  b  =  6.0  x  102,  b  =  6.0  x  101  and  3  —  1.8  X  10-5, 
3  =  1.8  x  10-4.  The  models  with  3  —  0  and  b  —  6.0  x  102  and  3  =  1.8  x  10~4. 
g(s)  =  0  show  considerable  increases  in  open-loop  damping  (especially  at  the  higher 


modes)  as  illustrated  in  Table  2. 


The  LQR  problem  for  the  approximating  system  (4.10)  is  solved  (by  Potters 
method)  and  approximating  functional  gains  A’(v(-),  A'3X  (•)  and  A'4V(-.-)  are 

constructed  by  the  standard  G'alerkin  scheme  (see  [M]  for  details).  Convergence  of 
these  sub-optimal  functional  gains  as  A’  — >  oc  to  the  optimal  gains  can  be  established 
for  3  >  0  and  g(s)  =  0  by  methods  similar  to  those  in  [G].  If  ^(s)  is  defined  by  (4.5), 
then  this  convergence  can  be  established  by  a  modification  of  Ito's  recent  results  [I]. 
Nothing  is  known  about  the  classical  thermo-elastic  problem.  We  treat  this  problem 
numerically. 

Example  4.1.  In  this  problem  we  set  <j{s)  =  0  —  3  and  compute  the  gains  AjN(-), 


( • )  and  A'3  (•).  Note  that  the  finite  dimensional  system  (*1.10)  is  exponentially 
stable  (all  open-loop  poles  are  in  the  left  half  plane)  so  that  these  sub-optimal  func¬ 
tional  gains  do  exist.  Figures  2,  3  and  *1  contain  the  plots  of  A';N(-),  i  =  1.2,3  for 
i V  =  24,  28,  32  and  36.  Although  these  plots  are  similar  to  those  found  in  [GRT], 
and  for  the  boundary  conditions  in  [GRT]  one  can  prove  convergence  of  the  gains, 
it  is  not  clear  from  Figures  3-4  that  these  gains  are  converging.  We  will  return  to 
this  issue  later. 

Example  4.2.  In  this  problem  we  consider  the  thermo- viscoelastic  problem  defined 
by  3  =  0  and  b  =  6.0  x  101,  b  —  6.0  x  1U\  Figures  5,  6  and  7  contain  the  plots  of 
A',A(-),  1  =  1.2,3  for  :Y  =  24,  28,  32  and  36  where  b  =  6.0  x  101.  Figures  8,  9  and 
10  contain  the  same  plots  for  6  =  6.0  x  102.  Observe  three  important  features:  i) 
For  b  =  6.0  x  101  the  system  is  more  heavily  damped  than  for  b  =  6.0  x  102:  ii)  the 
functional  gains  for  the  case  b  =  6.0  x  101  are  smoother  than  for  6  =  6.0  x  1 02 ;  and 
Hi)  the  convergence  of  the  functional  gains  is  faster  in  the  problem  with  the  most 
open-loop  damping.  The  plots  of  I\^{-,-)  for  b  =  6.0  x  101  and  N  =  24,  28,  32  and 
36  are  shown  in  Figures  11-14.  Note  that  these  gains  converge  rather  rapidly.  This 
was  typical  of  all  the  thermo- viscoelastic  runs. 

Example  4.3.  In  this  problem  we  consider  only  the  thermo-elastic  model  with 
Kelvin-Voigt  damping  so  that  <?(s)  =  0  and  3  —  1.8  x  10-4  and  3  =  1.8  x  10-5. 
Figures  15,  16  and  17  contain  the  plots  of  A',a(*).  i  =  1.2.3  for  A’  -  24.  28.  32  and  36 


•15 


for  3  =  1.8  x  10-4,  and  Figures  18,  19  and  20  show  the  same  plots  for  3  =  1.8  x  10~5. 
Observe  that  as  in  Example  -1.2  above,  the  damping  affects  the  smoothness  and  rate  of 
convergence  of  the  functional  gains.  Also  note  that  when  3  =  1.8  x  10-5  the  first  elas¬ 
tic  mode  has  the  same  damping  factor  (i.e.,  open-loop  pole  is  —  0.00039  +  6. 58i )  as  the 
problem  in  Example  4.2  with  b  —  6.0  x  102  (i.e.,  open-loop  pole  is  —0.00038  +  6.5Sz). 

Figures  21,  22  and  23  contain  the  plots  of  A  j,2(  • )  for  the  thermo-elastic  and  thermo¬ 
viscoelastic  models  and  the  thermo-clastic  model  with  Kelvin-Yoigt  damping  for  var¬ 
ious  values  of  b  and  3.  The  functional  gains  for  the  thermo- viscoelastic  problem 
"converge"  to  the  functional  gains  for  the  thermo-elastic  problem  as  b  — ♦  dc.  The 
same  convergence  applies  as  the  Kelvin- Voigt  parameter  3  —*  0.  As  illustrated  in 
Figure  23.  the  functional  gain  A'j^F)  computed  by  using  the  thermo-elastic  model  is 
the  same  as  the  gain  computed  by  using  the  same  model  with  Kelvin-Yoigt  damping 
parameter  3  =  1  .S  x  10”' .  Likewise,  if  b  ~  6.0  x  104  in  the  thermo-viscoelastic  model, 
then  the  functional  gain  A’j!2( • )  is  also  identical  to  the  gain  computed  by  using  the 
thermo-elastic  model.  This  “basic"  observation  applies  to  A',32(-)  for  i  —  1,2,3  and  for 
A',64(-)  also.  In  particular,  as  shown  in  Figure  24,  the  functional  gain  A'f4(-)  computed 
from  the  thermo-elastic  model  is  remarkably  close  to  the  functional  gain  computed  by 
adding  Kelvin-Yoigt  damping.  Although  convergence  of  the  functional  gains  for  the 
thermo-elastic  model  with  Kelvin-Yoigt  damping  is  assured  by  theory,  Figures  23  and 
24  illustrate  that  this  convergence  can  be  extremely  slow  if  the  damping  parameter 
3  is  small. 
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Figures  23  and  21  also  raise  the  possibility  that  even  in  the  thermo-elastic  problem 
(no  structural  damping)  where  there  is  no  proof  of  convergence,  the  gains  might 
converge  if  X  is  sufficiently  large.  The  jagged  nature  of  h\s  (■)  occurs  because  of  the 
small  damping.  Moreover,  comparing  the  open-loop  and  closed-loop  poles,  one  sees 
that  the  optimal  feedback  law  introduces  considerable  damping  at  all  frequencies. 

Figure  25  compares  the  first  S  open-loop  polos  (  +  )  to  the  closed-loop  poles  (*)  for 
the  thermo-elastic  model.  It  is  interesting  to  note  that  the  same  pattern  holds  for 
all  of  the  higher  poles  (see  Table  3  below),  and  the  thermal  damping  alone  does  not 
appear  to  be  a  major  aid  in  controlling  the  higher  modes. 


Mode 

Open-Loop 

Closed-Loop 

1 

-0.0000054  +  6.5Sz 

—  0. 107S901  +  6.5S? 

2 

— 0.0000060  +  13.17? 

-0.0552142  +  13. 17? 

3 

—  0.000007 1  +  19.79? 

—  0.079225 1  +  19.79? 

4 

-0.00000S6  +  26.47? 

-0.0953S6S  +  26.47? 

5 

-0.0000106  +  33.20? 

-0.029S920  +  33.20? 

6 

-0.0000129  +  40.02? 

-0.1097076  +  40.02? 

i 

-0.0000157  +  46.93? 

-0.0264249  +  46.93? 

S 

-0. 0000 1S9  +  53.95? 

—  0.09465S3  +  53.95? 

9 

-0.0000224  +  61.10? 

-0.0742450  +  61.10? 

10 

-0.0000263  +  6S.39? 

—0.0549855  +  68.39? 

15 

-0.0000491  +  107.5? 

-0.0707290  +  107.5? 

20 

-0.0000662  +  151.5? 

-0.0731 157  +  151.5? 

TABLE  3.  OPEN-LOOP  YS.  CLOSED-LOOP  POLES 


THERMO- ELASTIC  MODEL 
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As  illustrated  by  Figures  26  and  27,  the  first  S  closed-loop  poles  for  the  moderately 
damped  ( Kelvin- Voigt )  thermo-elastic  and  the  thermo- viscoelastic  models  are  essen¬ 
tially  the  same  as  the  corresponding  closed-loop  poles  for  the  thermo-elastic  model. 

In  view  of  these  numerical  results,  we  first  conjectured  that  although  no  theoretical 
results  exist  to  prove  the  existence  of  the  optimal  feedback  gain  for  the  classical 
thermo-elastic  model  (with  boundary  conditions  (4.3)),  existence  and  convergence  of 
the  suboptimal  gains  do  hold.  Since  the  damping  for  the  realistic  model  is  so  small, 
the  numerical  evidence  provided  by  the  previous  plots  is  not  strong.  In  order  to  test 
this  conjecture  we  investigated  the  convergence  of  the  gains  for  the  problem  with 
"artificial"  parameters.  In  particular,  we  considered  the  equations  in  dimensionalless 
form 

^2  cP  (} 

do(t,x)  =  ad-0(t.x)-c‘-tl^xy{t,X) 

with  b(x)  as  above  and  c2  —  7  =  1.0.  The  uncontrolled  system  has  considerable  damp¬ 
ing.  The  first  two  elastic  modes  have  eigenvalues  -.0939  +  3. lot  and  —.4066  +  6.47?', 
respectively.  The  functional  gains  for  this  set  of  parameters  are  slightly  smoother 
than  the  corresponding  gains  for  the  aluminum  rod.  Figures  28,  29  and  30  show  the 
“convergence"  of  these  gains  for  ,V  —  24,28,32  and  36.  It  appears  that  the  scheme 
may  produce  convergent  gains  (at  least  for  7\’jN  (•)  and  Kj  (•));  however,  the  behaviour 
of  h'2  ( ■ )  is  not  clear-cut.  Based  on  our  numerical  experience,  we  do  conjecture  that 
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the  classical  thermo-elastic  model  with  Dirichlet  boundary  conditions  is  in  fact  stabi- 
lizable.  However,  we  also  conjecture  that  the  computational  scheme  does  not  produce 
convergent  gains  (we  believe  that  this  scheme  does  not  preserve  stabilizability  uni¬ 
formly).  We  have  not  been  able  to  provide  proofs  of  these  conjectures,  and  we  leave 
them  as  open  problems. 
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